Information theoretic perspective on genome clustering.

Shannon's information theoretic perspective of communication helps one to understand the storage and processing of information in one-dimensional sequences. An information theoretic analysis of 937 available completely sequenced prokaryotic genomes and 238 eukaryotic chromosomes is presented. Information content (Id) values were used to cluster these chromosomes. Chargaff's second parity rule i.e compositional self-complementarity, an empirical fact is observed in all the genomes, except for the proteobacteria Candidatus Hodgkinia cicadicola. High information content, arising out of biased base composition in all the 14 chromosomes of Plasmodium falciparum is found among two other genomes of prokaryotes viz. Buchnera aphidicola str. Cc (Cinara cedri) and Candidatus Carsonella ruddii PV. Despite size and compositional variations, both prokaryotic and eukaryotic genomes do not deviate significantly from an equiprobable and random situation. Eukaryotic chromosomes of an organism tend to have similar informational restraints as seen when a simple distance based method is used to cluster them. In eukaryotes, in certain cases, Id values are also similar for the two arms (p and q arm) of the chromosomes. The results of this current study confirm that the information content can provide insights into the clustering of genomes and the evolution of messaging strategies of the genomes. An efficient and robust Perl CGI standalone tool is created based on this information theory algorithm for the analysis of the whole genomes and is made available at https://github.com/AlagurajVeluchamy/InformationTheory.

Understanding the virulence of the entero-aggregative E. coli O104:H4.

O104:H4 is a new strain of E. coli that has caused an outbreak in Germany. It was isolated from patients with bloody diarrhoea and Haemolytic Uremic Syndrome (HUS). BGI (www.genomics.cn) sequenced and assembled this new strain. It was reported to show resistance to a number of drugs that are toxic to other E. coli and causes serious complications during infections, which ultimately lead to death. Multi-drug resistance and high virulence of this strain is thought to be acquired from different sources, by horizontal gene transfer. A total of 38 prophage elements were detected from the new strain by using three computational tools viz., DRAD, Prophage Finder and PHAST. Analysis on these prophage elements shows a number of virulence proteins like Shiga toxin and multi-drug resistance protein encoding genes. The high virulence of the strain could be attributed by the prophage elements acquired from its micro environment.

Characterization of DLP12 Prophage Membrane Associated Protein: HolinGFP.

Lysis cassette genes from phages determine the final lytic event of the host cells. The lysis cassette genes are conserved in phages and prophages. The membrane associated holin from DLP12 prophage, available as a GFP fusion construct, was shown to be overexpressed, using confocal microscopy analysis, in bacterial cells. The protein expression caused cell death in E. coli AG1 strain suggesting the protein was functional. The His-tag HolinGFP protein was purified using cobalt affinity column and was eluted in the presence of different non-ionic detergents DDM (n-dodecyl-ß-d-maltoside), LDAO (Lauryldimethylamine-oxide), OG (n-octyl ß-d-glucopyranoside) and C12E9 (dodecyl nonaoxyethylene ether). HolinGFP existed predominantly as a dimer in LDAO in Superdex S200 gel filtration chromatography. Circular dichroism and fluorescence spectroscopy of the fluorescent HolinGFP in all four detergents (C12E9, DDM, LDAO, and OG) confirmed the folded state. Both dithiobis succinimidyl propionate and gluteraldehyde crosslinking revealed the existence of higher order oligomers and dimers. HolinGFP has been functionally and biophysically characterised and is being explored for crystallographic structure determination.

Molecular modeling and conformational analysis of native and refolded viral genome-linked protein of cardamom mosaic virus.

The viral genome-linked protein (VPg) of Potyviruses is covalently attached to the 5' end of the genomic RNA. Towards biophysical characterization, the VPg coding region of Cardamom mosaic virus (CdMV) was amplified from the cDNA and expressed in E. coli. Most of the expressed VPg aggregated as inclusion bodies that were solubilized with urea and refolded with L-arginine hydrochloride. The various forms of CdMV VPg (native, denatured and refolded) were purified and the conformational variations between these forms were observed with fluorescence spectroscopy. Native and refolded CdMV VPg showed unordered secondary structure in the circular dichroism (CD) spectrum. The model of CdMV VPg was built based on the crystal structure of phosphotriesterase (from Pseudomonas diminuta), which had the maximum sequence homology with VPg to identify the arrangement of conserved amino acids in the protein to study the functional diversity of VPg. This is the first report on the VPg of CdMV, which is classified as a new member of the Macluravirus genus of the Potyviridae family.

Identification of prophages in bacterial genomes by dinucleotide relative abundance difference.

BACKGROUND: Prophages are integrated viral forms in bacterial genomes that have been found to contribute to interstrain genetic variability. Many virulence-associated genes are reported to be prophage encoded. Present computational methods to detect prophages are either by identifying possible essential proteins such as integrases or by an extension of this technique, which involves identifying a region containing proteins similar to those occurring in prophages. These methods suffer due to the problem of low sequence similarity at the protein level, which suggests that a nucleotide based approach could be useful. METHODOLOGY: Earlier dinucleotide relative abundance (DRA) have been used to identify regions, which deviate from the neighborhood areas, in genomes. We have used the difference in the dinucleotide relative abundance (DRAD) between the bacterial and prophage DNA to aid location of DNA stretches that could be of prophage origin in bacterial genomes. Prophage sequences which deviate from bacterial regions in their dinucleotide frequencies are detected by scanning bacterial genome sequences. The method was validated using a subset of genomes with prophage data from literature reports. A web interface for prophage scan based on this method is available at http://bicmku.in:8082/prophagedb/dra.html. Two hundred bacterial genomes which do not have annotated prophages have been scanned for prophage regions using this method. CONCLUSIONS: The relative dinucleotide distribution difference helps detect prophage regions in genome sequences. The usefulness of this method is seen in the identification of 461 highly probable loci pertaining to prophages which have not been annotated so earlier. This work emphasizes the need to extend the efforts to detect and annotate prophage elements in genome sequences.

Subclassification and targeted characterization of prophage-encoded two-component cell lysis cassette.

Bacteriophage induced lysis of host bacterial cell is mediated by a two component cell lysis cassette comprised of holin and lysozyme. Prophages are integrated forms of bacteriophages in bacterial genomes providing a repertoire for bacterial evolution. Analysis using the prophage database (http://bicmku.in:8082) constructed by us showed 47 prophages were associated with putative two component cell lysis genes. These proteins cluster into four different subgroups. In this process, a putative holin (essd) and endolysin (ybcS), encoded by the defective lambdoid prophage DLP12 was found to be similar to two component cell lysis genes in functional bacteriophages like p21 and P1. The holin essd was found to have a characteristic dual start motif with two transmembrane regions and C-terminal charged residues as in class II holins. Expression of a fusion construct of essd in Escherichia coli showed slow growth. However, under appropriate conditions, this protein could be over expressed and purified for structure function studies. The second component of the cell lysis cassette, ybcS, was found to have an N-terminal SAR (Signal Arrest Release) transmembrane domain. The construct of ybcS has been over expressed in E.coli and the purified protein was functional, exhibiting lytic activity against E.coli and Salmonella typhi cell wall substrate. Such targeted sequence- structure-function characterization of proteins encoded by cryptic prophages will help understand the contribution of prophage proteins to bacterial evolution.

A simulation model of Escherichia coli osmoregulatory switch using E-CELL system.

BACKGROUND: Bacterial signal transduction mechanism referred to as a "two component regulatory systems" contributes to the overall adaptability of the bacteria by regulating the gene expression. Osmoregulation is one of the well-studied two component regulatory systems comprising of the sensor, EnvZ and the cognate response regulator, OmpR, which together control the expression of OmpC and OmpF porins in response to the osmolyte concentration. RESULTS: A quantitative model of the osmoregulatory switch operative in Escherichia coli was constructed by integrating the enzyme rate equations using E-CELL system. Using the substance reactor logic of the E-CELL system, a total of 28 reactions were defined from the injection of osmolyte till the regulated expression of porins by employing the experimental kinetic constants as reported in literature. In the case of low osmolarity, steady state production of OmpF and repression of OmpC was significant. In this model we show that the steady state - production of OmpF is dramatically reduced in the high osmolarity medium. The rate of OmpC production increased after sucrose addition, which is comparable with literature results. The relative porin production seems to be unaltered with changes in cell volume changes, ATP, EnvZ and OmpR at low and high osmolarity conditions. But the reach of saturation was rapid at high and low osmolarity with altered levels of the above components. CONCLUSIONS: The E-CELL system allows us to perform virtual experiments on the bacterial osmoregulation model. This model does not take into account interaction with other networks in the cell. It suggests that the regulation of OmpF and OmpC is a direct consequence of the level of OmpRP in the cell and is dependent on the way in which OmpRP interacts with ompF and ompC regulatory regions. The preliminary simulation experiment indicates that both reaching steady state expression and saturation is delayed in the case of OmpC compared to OmpF. Experimental analysis will help improve the model. The model captures the basic features of the generally accepted view of EnvZ-OmpR signaling and is a reasonable starting point for building sophisticated models and explaining quantitative features of the system.